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ABSTRACT 

The study of young stellar populations has revealed that most stars are in binary 
or higher order multiple systems. In this study the influence on the stellar initial 
mass function (IMF) of large quantities of unresolved multiple massive stars is investi- 
gated by taking into account stellar evolution and photometrically determined system 
masses. The models where initial masses are derived from the luminosity and colour 
of unresolved multiple systems show that even under extreme circumstances (100% 
binaries or higher order multiples) the difference between the power-law index of the 
mass function of all stars and the observed mass function is small (^ 0.1). Thus, if 
the observed IMF has the Salpeter index a = 2.35 then the true stellar IMF has an 
index not flatter than a = 2.25. Additionally, unresolved multiple systems may hide 
between 15 and 60% of the underlying true mass of a star cluster. While already a 
known result, it is important to point out that the presence of a large number of un- 
resolved binaries amongst pre-main-sequence (PMS) stars induces a significant spread 
in the measured ages of these stars even if there is none. Also, lower-mass stars in a 
single-age binary-rich cluster appear older than the massive stars by about 0.6 Myr. 

Key words: binaries: close - binaries: general - stars: early-type - stars: evolution 
- stars: formation - stars: luminosity function, mass function 



1 INTRODUCTION 



Most late-type stars ( i$ IMq) are born in clusters of binary systems with companion masses distributed randomly from the stel- 
lar in itial mass function (IMF; Goodwin fc Kroupal2005l : Goodwin et aLlBoOT : Goodwin fc Whitworthl2007 : Kouwenhoven et al 
20081 ). Subsequent dynamical encounters i n the cluster evolve the stellar population to that observed in the Galactic field 



IPuguennoy fc Mayor 1991 : Kroupa 1995a ^QE; Kroupa fc Bouvier 20031 ). If massive (^fewM©) stars were also to be born 



with companions distributed randomly from the IMF, then the vast majority of O stars would have M-dwa rf companions while 
mass ive equal-mass binaries would be extremely unlikely. For random sampling from a canonical IMF ( Weidner fc Kroupal 



200 



3) for any O stai0it would be expected that 0.2% of the companions are O stars (m ^ IGM©), 2.2% B stars (2.7Mq 16Mq) 



2.1% A stars {1.7Mq <m<2.7MQ), 3.7% F stars (I.IM© < m < 1.7AfQ), 3.4% G stars (0.85Mq < m < I.IMq), 12.0% K stars 



* E-mail: Carsten.Weidner@st-andrews.ac.uk 
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^ The percentages are valid for any given star but if the "companion" is more massive than the primary they would switch places and 
the ratios would change. If we constrain ourselves to O stars as primaries it does not matter if the companion is more massive than the 
primary because they would be both O stars. 
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(O.SMg 



,<0.85Mq) and 76.3% M stars (0.08 Mg 



,'^0.5Mq). These M-dwarfs (as well as the K and at least the G 



stars, altogether more than 90% of all companions) would b e more or less invisible. However, the observed fraction of O stars in 



mass i ve multiple systems l ies between at least 20 and 80 % ( Garmanv et al, 
20061 : iKiminki et aLllioO?! : iLucvlbood : Upai et al.ll2007l : ISana et allboOSl : 



198C 



Garcia fc M ermilliod '2001 
Kouwenhoven et al.i,2007i : .Turner et al 



De Becker et al 



2003) and 



thus strongly contradicts secondaries which are randomly taken from the IMF, at least for massive primaries. Unfortunately, 
the mass ratio {q = !^££22;is2L) distributioiH for massive stars is observationally not well explored for several reasons: Massive 
stars have very broad spectral lines, making radial velocity studies a challenge. Even rather massive and close companions 
are easy to miss. And young stellar nurseries are usually deeply embedded in gas and dust. Furthermore, massive stars are 
predominately found in massive clusters. So the issue of crowding and fast dynamical evolution arises and therefore it is not 
clear if the gr-distribution is primordial. As massive stars either sink to the centre of star cl usters through energy equipartition 



GrebellboOSi: 



Stolte et al 



Bonnell et al. 



1998 



2003 



20061 : iFleck et allbood : IVesperini et allbood : Izhao et allbood ) or are (exclusively?) formed there 



2004 ) ■ dynamical processes in the cl uster centre can expel (massive) stars from the cluster and/o r 



exchange binary partners, thereby changing the g-distribution ( Pflamm- Altenburg fc Kroupa 20061 : Pfalzner fc Olczak 2007 ). 
This picture is strongly s upported by the fact that "field" and ru naway O stars have a far lower binary fraction than O stars in 
clusters and associations ( Turner et al. 20081 : Mason et al. 20081 ). A recent study of triple and quadruple system by TokovininI 
1 20081 ) might also be used as an argument against random pairing of stars as the author finds a surprisingly large amount of 
these syste ms in which all sta rs are of similar mass. Whether or not these results are applicable to this study is a matter of 
debate. The Tokovinin ( 20081 ) sample consists to about 80% of primary stars below 5 Mq and it is not restricted to young 
objects. If the author's claim is correct that properties like orbital periods and eccentricities of old multiple systems are not 
consistent with A^-body calculations his study would be an important argument against random pairing of multiple stars. It 
should be noted here that in a cluster of the richness of the Orion Neb ula cluster stellar dynamics may change the mass func- 
tion and binary properties in the core significantly in less than 10^ yr ( Pflamm- Altenburg fc Kroupal 200d : Pfalzner fc Olczak 
20071 ). The study of the q-distribution in the present paper is of stars at the stage when they are observed and may thus be 
representative of dynamically evolved distributions. 



The influ ence of a substantial fraction of unresolved binary systems on the lo w-mass stellar mass function (MF) has been 
wefl e xplor ed (jSaear fc Richtlejll99ll : iKroupa et al]ll99ll . Il993l : [Kroupa" 1995a''d'; 'Malkov fc Zinne ckerllioOll : iThies fc Kroupal 
2007 . 20081 ) . Beginning from the lowest mass stars and brown dwarfs. Thies fc Kroupa (2007.. .20081 ) corrected for un-resolved 



companions in such system s and showed that the IMF is discontinuous near 0.1 Mq. In an extensive study of late-type stars, 
Kroupa et al.l (119911 . Il993l ) showed that the mass function of these stars is signiflcantly affected by unresolved systems. They 
used this result to formulate the standard or canonical-IMF with little or no evidence for systematic variations with star 
forming conditions (|Kroupall200ll . for more details on the canonical IMF: see App endix [A|l . Interestingly, there seems now 
to be strong evidence for a universal companion mass function for solar type stars ( Metchev fc Hillenbrandl 20081 ). But only 
a ver y limited number of stu dies have addressed how the large proportion of binaries affects IMF derivations for massive 
stars. ISagar fc Richtlerl l|l99ll ) studied the MFs of five young star clusters in the Large Magellanic Cloud (LMC) including 
the influence of optical and physical binaries on the derived slope of the MFs in the range of 2-14 Mq by Monte Carlo 
simulations of cluster and field populations with random pairing over this mass int e rval. They found a very strong steepening 
of the observed MF in the case of a binary fraction larger than 50%. IVanbeverer] (|l982l l investigated the IMFs of primaries 
and single stars analytically up to stars above 100 Mq but did not explicit l y disc uss any possible effects of unresolved binaries 
on the observed IMF. Only rather recently, the study by Mai'z Apellaniz ( 20081 ) explores the biases on the IMF introduced 
by unresolved binaries for massive stars. Ther e a different col o ur is used {U — V) while in this contribution we resort to 
B — V, and only random pairing is studied bv lMai'z Apellaniz (12008', see § 12.11 for details on binary pairing) while here we 
consider different viable q-distributions. This contribution can therefore be seen as an extension of thes e previous studies 
( Kroupa et al]ll99ll . 1 19931 : ISagar fc Richtleijll99ll: [ Malkov fc Zinneckej [20011 : iThies fc Kroupal l2007l . I2OO8I ) to higher masses 



as well as a comparison to the lMaiz Apellaniz! (|2008l ) study. 



2 THE MODEL 

In what follows, dN = ^(m) dm is the number of stars in the mass interval m to m 4- dm. 

For all methods the IMF, £,{m), described in Appendix]^ is used as the input IMF. For every model the percentage of 
binaries, triples and quadruples is assigned as an initial condition and then one of the following pairing algorithms is used to 
obtain the masses of the stars. The primary, secondary etc. is always the most-massive, second massive, etc., companion. 



The q-distribution for r andom pairing is roughly flat, showing little preference for any particular q-value jKroupal Il995d : 
iGrether fc Lineweave3l2006l) . 
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2.1 Pairing algorithms 

As the initial properties of (massive) binaries and liiglier-order multiple systems are generally unknown, several difTerent 
methods of pairing stars were used. 

• Random pairing (RP): A given number of stars is taken fully randomly from the input IMF over the mass range of 0.08 
to 150 Mq and randomly assigned to pairs or higher-order multiple systems. 

• Special pairing (SpP): Again a given number of stars is taken randomly from the input IMF over the mass range of 0.08 
to 150 Mq and randomly grouped together. But in the cases where the primary exceeds 2 Mq the secondary has to fulfil 
a minimum mass criterion as described in Appendix IBll If the secondary fails this test it is exchanged with another star 
from the random list until the minimum mass is reached or exceeded. In the case that the new secondary is more massive 
than the primary the two exchange places and the secondary becomes the new primary and the old one the new secondary. 
This new secondary still has to fulfil the minimum mass criteria or it is exchanged with another star. It is possible that no 
such secondary is available in the list anymore. Then the minimum mass is lowered by 20% and if this is still insufficient, 
the star is kept as a single. In the numerical experiments presented here no such case was encountered. In the case of triple 
or quadruple systems, if the secondary component is also larger than 2 Mq the tertiary component is searched for with the 
same criterion as the secondary before (ie the secondary is treated as the "primary" in the secondary-tertiary sub-system). 
The same procedure is extended to the quartiary companions (ie. the quartiary becomes the "secondary" of the tertiary) . 

• System pairing (SyP): A given number of system masses is randomly taken from an input (not the canonical) IMF with 
0.16 ^ rrisys ^ 300 Mq and with a change of slope from 1.3 to 2.35 at 1 Mq instead of 0.5 Mq. The lower and upper mass 
limit for stars is always kept at 0.08 and 150 Mq. If the system mass is below 10 Mq the system is randomly split into two, 
three or four components assuming a uniform mass-ratio distribution being equivalent to random pairing on this mass interval. 
But for systems with a mass that exceeds 10 AIq the g-values used to split the systems into components are biased towards 
high-g- values as described in Appendix lB2l This model is motivated by a possib le merger-origin of massive stars but could also 



nign-g- values as described m Appendix lczl inis model is motivated by a possip le merger-origin oi massive stars but could also 
represent the outcome of the competitive accretion model i Bonnell et al.lToQSl ). In order to have a smooth transition between 



the random mass ratios and the preferential formation of massive companions for the SyP model the change is necessary at 
10 Mq and not at 2 Mq as in the case of the SpP model. 

All pairing mechanisms preserve a lower mass limit of 0.08 Mq for all components, because stars and brown dwarfs rarely 



pair up dThies fc KroupalboO?! . I2OO8I ) 



After the stars are assigned into their systems, the mass function for all stars, for the primaries, secondaries, tertiaries (if 
applicable) and quartiaries (if applicable) is calculated along with the IMF for the systems. The method used to measure the 
slopes is presented in § 12.31 The direct results of these calculations are presented in § 13.11 while the calculations additionally 
applying stellar evolution (as described in § 12. 2p are discussed in § 13.21 



2.2 Stellar evolution 



Stellar observations usually do not result directly in stellar masses and ages but magnitudes and colours. Therefore, the stars 
produced by the pairing mechanisms described in Section [2TT] are assigned a random age with a Gaussian distribution between 
0.01 and 4.5 Myr with a median of 2 Myr and a dispersion of 0.5 Myr. Stars in the same multiple system are assumed to be 
coeval. Using a large grid of stellar model^ this age is then used to evolve the stars and obtain the evolved masses, effective 
temperatures and r adii. These values are further used to calculate the Johnson V and B band magnitudes with the help of 
Kurucj (|l992l ') and [Bergeron et al. l|l995l ) bolometric corrections using a program kindly provided by J. Hurley (see iHurlevI 
2OO3I ). The luminosities of the individual stars are then combined to a single V magnitude and B ~V colour of the unresolved 
multiple and searched for in a large table of V and B ~ V values made with the same models but for truly single stars and 
with an age range between 0.01 and 100 Myfl In this manner so-called pseudo masses and pseudo ages for the unresolved 
multiples are derived as if they were single stars. Note that the derived pseudo masses are the zero-age masses. 

Note that photometric mass and age determination is know to be inaccurate in comparison to spectroscopic measurements 
i MassevHiooi ). But as spectra are far more expensive in terms of observational time the photometric method is still widely 
in use. Furthermore, for OB stars a spectroscopic analyses does not automatically solve the problem of un-resolved multiples 
as O stars have very broad lines which easily hide even massive compa nions. Only with elaborate multi-epo ch spectroscopic 
observations solid constrains can be put on the multiplicity of O stars ( Sana et al. 20081 : Schnurr et al. 20081 '). 



3 iMevnet fc Maedej | |2003^ models for mas s ive stars. iHur lev et aL I ll2000h for the evolution after the pre-main sequence (PMS ) for 
stars below 9 Mq and iBurrows et al.l lll993h . iBurrows et al.l l|l997l 'l. IP'Antona fc Mazzitellil l|l997l ). IP'Antona fc MazziteTIi 1 1996 ) and 
iBehrend fc MaedeJ 1I2OOII ') for PMS evolution. For further details on the used models for stellar evolution see lWeidner et al.l 1 12008) . 
^ A larger age than the input age of 4.5 Myr is used to additionally show the effect of unresolved multiples on the age determination of 
stars. 
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Table 1. Standard deviation and bias for the estimation method of the exponent, for different sample sizes and output exponents. 
Standard deviation 



# stars 


a = 1.5 




a = 2.3 




a = 5.0 




10^ 


0.026 




0.043 




0.127 




10-* 


0.009 




0.014 




0.042 




10^ 


0.003 




0.004 




0.012 




Bias 














#: stars 


a = 1.5 




a = 2.3 




a = 5.0 




10^ 


3 X 10" 


4 


3 X 10- 


3 


1 X 10- 


2 


10* 


4 X 10" 


5 


3 X 10- 


4 


2 X 10" 


3 


10^ 


1 X 10- 


5 


3 X 10- 


5 


1 X 10" 


4 



2.3 A bias- free fitting method 



The common method to analyse power-law distributed data is to group them in bins of constant logarithmic size and then 
perform a linear regression to measure the power-law index a. However, this method c an lea d to strongly biased results 
i Mafz Apellaniz fc Ubeda 2005 ). Therefore here the method of Maschberger fc Kroupa ( 20081 ') is used. It is based on the 

estimates for a obtained with the maximum likelihood method. 

The maximum likelihood estimate of the power-law exponent (see Aban et al. 2006h is calculated by minimising the 
logarithmic likelihood function of the data {Xi}"^i,where Xi = rrii is the mass of star i. The estimated exponent uml is then 
the solution of 

Z^-^ML log z - yi-^Afi log Y 



+ n- 



T 



0, 



with the smallest data point, Y — minX^, the largest data point, Z 
T — logXi. The bias-free estimate follows with 

a = -(cvAfL - 1) + 1. 



(1) 



maxXj, and the sum of the logarithms of all data, 



(2) 



The standard deviation and bias of the method were investigated by using a Monte-Carlo sample of 1000 synthetic data 
sets, generated with a lower limit of 1.3 Mq, an upper limit of 150 Mq, different sizes of the data set and different exponents. 
For each data set the exponent was estimated and from the sample of estimates the standard deviation and the bias (difference 
between the average of the estimated exponents and the input value) were calculated. The results are summarised in Table 
[1] The bias is in all cases negligible. The standard deviation decreases with increasing size of the data set and increases with 
increasing exponent. For the sample sizes used in this work a conservative estimate of the standard deviation is 0.05. 



3 RESULTS 

In § 13.11 the results of simulating binaries, triples and quadruples are described under the assumption that all stellar masses 
are known perfectly. In i; I3.2l the stars are evolved with the use of stellar evolution models and the combined magnitudes and 
colours are used to model the process of observing stellar populations containing unresolved multiple systems. 



3.1 Models without stellar evolution 

Shown in Tab. [2] are the parameters and results for a series of models. A number of abbreviations are used in the Table. 
"100% bin" marks models made of 100% binaries, "100% trip" are 100% triples and "100% quad" are 100% quadruples. The 
corresponding IMF slopes are: QaiistarsEI which accounts for all stars individually (ie. the true single-star IMF), asystcms which is 
the mass function for the combined masses of all multiple systems, aprimaries which is the IMF of the most massive component 
of a system, asecondarios for the second most massive component of a system, atertiarios for the third massive object (only in 
triples and quadruplets) and finally aquartiarics for the fourth component in a system if there is any (only quadruplets). 

All models are calculated without stellar evolution, meaning they are treated as if all masses are known precisely. For 
each model 10^ stars are generated, each model taking only about a few minutes of computational time on a standard desktop 
PC. 



^ All IMF slopes in this work are obtained by a bias-free fitting method discussed in § 12.31 for all stars above 1 Mq . In the case of strong 
non-linear features in the MPs the line was fitted to the stars above such features. The standard deviation in a is of the order of 0.05 
dex. 
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Table 2. Properties of the models without stellar evolution. The pairing methods are described in Section l2.ll All stars are in 100% 
binaries, triples or quadruples, respectively. The slopes are measured by the bias-free method (see ^ 12.31 1 above 1.0 Mq unless noted 
otherwise. For each model a total of 10^ stars are generated. Column #1 denotes the model name, #2 the pairing method, #3 if the 
systems are binaries, triples or quadruples, #4 shows the slope for all stars, #5 the slope of the systems, #6 the slope of the primaries, 
#7 the slope of the secondaries, #8 the slope of the tertiaries and #9 the slope of the quartiaries. 



Model 


pairing method 


multiple 


^allstars 


^systems 


'^primaries 


secondaries 


'^tertiaries 


•^quartiaries 


1 


RP 


bin 


2.35 


2.40^ 


2.32 


3.71 






2 


RP 


trip 


2.35 


2.47I 


2.28 


3.64 


5.14 




3 


RP 


quad 


2.35 


2.54I 


2.26 


3.60 


4.96 


6.06 


4 


SpP 


bin 


2.35 


2.34I 


2.36 


2.042 






5 


SpP 


trip 


2.35 


2.38^ 


2.34 


2.13^ 


1.93 




6 


SpP 


quad 


2.35 


2.42^ 


2.31 


2.22^ 


1.78^ 


1.62^ 


7 


SyP 


bin 


2.35 


2.35 


2.60^ 


2.27 






8 


SyP 


trip 


2.35 


2.35 


2.57I 


2.26 


2.25 




9 


SyP 


quad 


2.35 


2.35 


2.55I 


2.25 


2.25 


2.25 


10 


q=l 


bin 


2.35^ 


2.35^ 


2.32 


2.36^ 







^ slope measured above 10.0 Mq, ^ slope measured above 2.0 Mq 
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Figure 1. Stylised flow chart how stellar evolution and unresolved binaries influence the observed mass function. 



As is seen in Tab. [2] and the first three panels of Fig. ICll the models using random sampling (RP Model 1 to 3) produce 
system IMFs which are somewhat steeper than the input IMF, and very steep IMFs for the companions, while the primary 
IMFs are slightly flatter. This can be understood as due to random pairing producing systems in which most of the massive 
stars have low-mass companions. These do not change the mass of the system significantly. But the low-mass systems tend to 
have partners of very similar mass, therefore additionally depleting the total number of low-mass systems and shifting them 
to higher masses. 

The special pairing mechanism explored here (SpP Model 4 to 6) has a different effect. The system IMF is also steepened 
but only for the case of 100% quadruples (SpP Model 6) and the companion star IMFs are flattened below the Salpeter-value 
for the secondaries, tertiaries and quartiaries for all three models. But for the primaries the IMF becomes slightly steeper. In 
the panels 4, 5 and 6 of Figs. ICll and IC2I this can be seen. This is expected since the SpP mechanism searches for massive 
companions for massive stars and therefore depletes the reservoir of massive stars remaining as primaries, and in doing so 
channels large amounts of massive stars into the companion star IMFs. 

A similar, but even more drastic result is obtained by system pairing (SyP Models 7 to 9, panels 7 to 9 in Figs. IC2l and 
C3|l . Here the primary IMFs are steepened and the secondary IMFs flattened in order to retain constant slopes for the system 
IMFs under the condition of having increasingly higher g- values for more massive systems. 

The rather flat slopes of the companion mass functions for SpP and SyP would be in agreement with the results of 
Tokovinin |2008l ) that for a large fraction of triple and quadruple systems the components are of rather similar mass. 



For reasons of comparison also an example of 100% equal-mass binaries for massive stars (for primaries > 2 Mq) is shown 
(Model 10 in Tab. [2] and panel 10 in Fig.[C3]|. Here in principle aU IMFs are parallel but again discontinuous features appear 
in the IMFs at the borders of the random pairing for the low-mass stars and the q = 1 pairing for the more massive ones. 



3.2 Models with stellar evolution 

In a second set of models the influence of the conversion of luminosities into masses is studied. The systems have ages assigned 
from a Gaussian distribution with a mean of 2 Myr and a standard deviation of 0.5 Myr and a cutoff at 0.01 and 4.5 Myr. 
With this age the stars are evolved according to the models described in Section [2. 2 1 The resulting stellar parameters are used 
to calculate luminosities and colours for each single star. These are merged to form an unresolved multiple system ensuring 
coevality of all companions. The resulting magnitudes and colours are searched for in a large table of colours/magnitudes of 
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Table 3. Like Tab.[2]but for the models with stellar evolution. For each model a total of 1 million stars are generated. Column #1 (model) 
denotes the model name, #2 (pairing method) the pairing method, #3 (multiple) if the systems are binaries, triples or quadruples, #4 
(c^allstars) shows the IMF slope for all stars above 1 Mq, #5 (cfsystoms) the slope of the systems, #6 (oprim) the slope of the primaries, 
#7 (osec) the slope of the secondaries, #8 (atort) the slope of the tertiaries, #9 (oquart) the slope of the quartiaries, #10 (Qobscrvod) 
the slope of the recovered "observed" stars, #11 (age fit) marks the percentage for how many of the unresolved stars the pseudo age 
of the system is recovered within 25% of the original age, #12 (total mass recovered) shows the percentage of the input mass in stars 
which was recovered in "observed" stars and #12 (age fit mass) is the recovered mass in stars for which the pseudo ages are within 25% 
of the real ages. The differences in the values for the slopes of the primaries, secondaries, tertiaries and quartiaries in comparison with 
Tab. [2] are due to the lower number of stars used for the models with stellar evolution, as these lower numbers significantly increase the 
error of the slope determination. 



Model 


pairing 


multiple 


'-^^allstars 


Q^systcms 






Otcrt 


C^quart 


'^observed 


age fit 


total mass 


age fit 




method 




















recovered 


mass 


Fit 0^ 


RP 


bin 


2.35 


2.40^ 


2.32 


3.71 






2.362 


50% 


84% 


59% 


Fit 1 


RP 


bin 


2.35 


2.40^ 


2.32 


3.70 






2.35^ 


50% 


85% 


43% 


Fit 2 


RP 


trip 


2.35 


2.47^ 


2.29 


3.63 


5.11 




2.343 


34% 


73% 


43% 


Fit 3 


RP 


quad 


2.35 


2.542 


2.25 


3.62 


4.98 


6.08 


2.333 


25% 


65% 


33% 


Fit 4 


SpP 


bin 


2.35 


2.36^ 


2.37 


2.07^ 






2.433 


49% 


80% 


44% 


Fit 5 


SpP 


trip 


2.35 


2.36^ 


2.34 


2.15^ 


1.98 




2.423 


32% 


68% 


30% 


Fit 6 


SpP 


quad 


2.35 


2.40^ 


2.32 


2.22^ 


1.83=* 


1.673 


2.433 


23% 


60% 


23% 


Fit 7 


SyP 


bin 


2.35 


2.35 


2.60^ 


2.26 






2.4l2 


43% 


70% 


25% 


Fit 8 


SyP 


trip 


2.35 


2.35 


2.56^ 


2.25 


2.24 




2.462 


22% 


53% 


10% 


Fit 9 


SyP 


quad 


2.35 


2.35 


2.542 


2.25 


2.25 


2.23 


2.432 


19% 


43% 


5% 


Fit 10 


q=l 


bin 


2.35^ 


2.352 


2.31 


2.372 






2.3l2 


45% 


68% 


10% 



^ This model is included for comparison only. It is identical to Fit 1 but all stars have a fixed input age of 2 Myr instead of the 

Gaussian distributed ages for the other models. 
2 slope measured above 10.0 Mq, 3 slope measured above 2.0 Mq 



single stars produced with the use of the same models as in Section [2.21 in order to find a pseudo zero-age mass and a pseudo 
age for the unresolved multiple as if it were a single star. The table covers a mass range between 0.01 Mq and 300 Mq with 
ages between zero and 100 Myr. In Fig. [T] the steps described above are visualised. As this procedure is rather CPU time 
consuming only 1 million stars are generated for each model in contrast to the 10 million in Section [3.11 Each model needs 
about 27 hours of CPU time on a state-of-the-art PC. 

The results and parameters for this second series of models are listed in Tab. |3] Additionally to the columns shown 
in Tab. [2] aobsorvod is added: the IMF slope of the mass function for all unresolved systems above 1 M© recovered by the 
simulated observations. In some models strong non-linear features occur in the output close to 1 Mq. In such cases a higher 
lower limit for the slope determination is used as indicated in the Table. The slope is determined in the same way as in § 13.11 
Furthermore, included is the column "age fit" in Tab. [3] which is the percentage of how many pseudo stars are recovered with 
a pseudo age within 25% of the real age of the system. Finally, the last two columns contain the recovered total mass fraction 
and the recovered mass fraction within the 25% age limit. In the first case the total mass of all recovered pseudo stars for 
one model is divided by the true stellar mass within the model and in the second case only the pseudo mass of the recovered 
pseudo stars with pseudo ages within 25% of their real ages is divided by the total true mass. Furthermore, for reasons of 
comparison. Model Fit is only included in Tab.[3]but not otherwise used. This model is similar to Model Fit 1 in all aspects 
but the age. Instead of an age spread all stars have the same input age of 2 Myr. 

Because of the fewer stars used here in comparison with i; 13. li the errors in the slope determination are larger than in the 
previous experiments without stellar evolution. 

In the first panel of Fig. IDll (Model Fit 1, random pairing, 100% binaries, see also the exemplary HR-Diagrams for this 
model in Fig. [2]) the input, primary, secondary and system IMF are the same as for Model 1 from Section [3. II but with a larger 
scatter because of the smaller number of generated stars. The derived observed IMF does not follow the (slightly) steeper 
system IMF but follows instead the input IMF of all stars. The steeping effect seen in the system IMF in comparison to the 
IMF of all stars is thus fully negated. This is due to the fact that while in the model without stellar evolution all masses are 
100% accounted for, in this model with evolution the low-mass companions of massive stars are "lost" because the primary 
dominates the light and only its mass is recovered. In the low-mass regime all stars are of rather similar mass and therefore 
the unresolved binaries change the luminosity of the system but not the {B — V) colour. This makes them appear as younger 
PMS stars, usually rather close to the mass of the primary. Both model series with mass- dependent mass-ratios (SpP and SyP, 
panels 4 to 9 in Figs. ID2l and lD3] ) have slightly steeper observed IMFs with a slope of about 2.45. But such a small deviation 
from the Salpeter-slope is within the error margin of measured high-mass IMF slopes. The final Model Fit 10 in which all 
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stars above 2 Mq reside in equal mass binaries influences the slope of the primary IMF and the observed IMF slightly. But 
such small differences (Aa = 0.04) are still within the errors of the here used slope determinations. 

In conclusion, the analysis presented here quite clearly shows that the observed IMF for massive stars remains mostly 
indistinguishable from the underlying true stellar IMF, which may be flatter by at most Aq = 0.1 even if all massive stars reside 
in quadruple systems. This would thus indicate that measured IMF indices are quite robust and that the Salpeter/Massey 
slope, a = 2.35, is not significantly affected by unresolved multiple stars above a few Mq. 
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Figure 3. HR-Diagrams for the systems and the observed pseudo stars for 17000 stars from model Fit 9. The solid line is a 2 Myr 
isochrone while the dotted line is a zero-age main-sequence. 



Table 4. Percentages of recovered ages in different age bins for a number of models. 



Model 


< 0.5 


0.5 - 1.5 


1.5 - 5.0 


5.0 - 10.0 


10+ 




Myr 


Myr 


Myr 


Myr 


Myr 


random ages (like Fit 1 to Fit 10), 












50% bin RP 


1.5 


19.7 


74.1 


1.6 


3.1 


0.5 Myr, 100% bin RP 


58.5 


22.1 


11.5 


3.5 


4.4 


1.0 Myr, 100% bin RP 


8.0 


65.2 


10.3 


9.0 


7.5 


2.0 Myr, 100% bin RP 


4.8 


28.9 


53.0 


4.3 


9.0 


3.0 Myr, 100% bin RP 


3.0 


10.3 


75.0 


4.8 


7.0 


Fit 1 


4.4 


31.8 


50.0 


4.6 


9.1 


Fit 4 


4.4 


31.8 


49.3 


4.9 


9.5 


Fit 7 


5.3 


38.0 


41.9 


3.8 


11.1 


Fit 10 


4.7 


33.9 


45.0 


4.7 


11.5 



3.2.1 Age spreads 

Interestingly, as is seen in the ninth column of Tab. |3] only about 40 to 50% of the pseudo stars are recovered with the right 
pseudo age (with an error of less than 25% compared to the true age). Actually, the discrepancy between assigned age and 
fitted age can be rather large. Therefore, large numbers of binaries (regardless of their masses) can mimic a wide spread of 
formation ages. In order to further explore this interesting result a series of models is calculated with 5 • 10^ binaries (0.08 to 
150 Mq for the stars), random pairing and the same input age. The distribution of the resulting recovered ages (using the 
same method as before in § I3.2|l for stars of all masses can be seen in Figs. |4]to [T] Fig. |4] features an input age of 0.5 Myr, 
Fig.[5]l Myr, Fig.|6]2 Myr and Fig. [7] 3 Myr. 

Using such a single input age (e.g. 1 Myr) as in Fig. [5] the distribution of the recovered output ages is severely distorted 
(dotted histogram in Fig. [5]). In Tab. |4] are shown the percentages of recovered ages in several age bins. Besides the models 
with single input ages also the Models Fit 1, Fit 4, Fit 7 and Fit 10 are included as well as a model with only 50% binaries 
and the same age range as the Models Fit 1 to Fit 10. 

In a real star cluster it is unlikely that all stars are truly of the same age. Therefore in Fig. [8] a distribution of input ages 
is used. Fig. [8]shows as an example the age distribution for Model Fit 1. While the solid histogram is the Gaussian input age 
distribution (distorted because of the logarithmic ordinate), the dotted histogram is the recovered age distribution. The mean 
of the recovered ages of all stars is slightly low er t han the true mean age (1.6 M yr instead of 2 Myr). This is consistent with 
the findings by Brandner fc Zinnecker (1997) and Preibisch fc Zinnecker 1 19991 ) that un-resolved binaries let a population 
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logio(age/Myr) 

Figure 4. Logarithm of the number of stars versus logarithm of the age for one million stars in 100% unresolved RP binaries (500000 
systems). Solid histogram: The input age of 0.5 Myr. Dotted histogram: The recovered output ages for stars of all masses. 




logio(age/Myr) 

Figure 5. Like Fig. |4]but with an input age of 1 Myr. 



appear younger than it is. Additionally, the wings of the recovered distribution are much wider - suggesting the presence of 
stars even up to an age of 100 Myr, despite the maximal inpu t age of about 4.5 Myr. Variability of thes e PMS stars and 
photometric errors make the age derivation even more difficult (IBriceno et al. 2007 : Hillenbrand et al. 20081 ). 

In Fig.|9]the same recovered output ages as in Fig.[8]are shown but split into two, for stars above 1 Mq {dotted histogram) 
and below {long-dashed histogram). It can be seen that the tail towards old ages is prominent in both populations but the tail 
on the left hand side, towards young ages, is more dominant for low-mass stars. However, the actual peak of the high-mass 
distribution is moved from 2 Myr to 1 Myr while the low-mass stars peak in between at about 1.6 Myr. The bulk low-mass 
stellar population therefore appears older by 0.6 Myr than the massive stars. 
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-10 12 
logio(age/Myr) 

Figure 6. Like Fig. |4]but with an input age of 2 Myr. 




logio(age/Myr) 

Figure 7. Like Fig. |4]but with an input age of 3 Myr. 



So, it is important to note that unresolved binaries easily lead to the impression of extended star-formation in a cluster and 
to the low-mass stars appearing older. Therefore, any age spread encountered in young stellar clusters needs to be interpreted 
very carefully, as there are al so other processes which may masquerade an age spread , e.g. older stars captured during the 



2006 



IPflamm-Altenburg fc Kroupall2o"o7l ) and/or unresolved multiple stars as 



formation of the star cluster (jFellhauer et al 
shown here. 

Figs. [10] and [n] show some details how the difference in the age determination comes about. The solid line in both 
figures refers to the evolution of an unresolved binary composed of equal age 1 Mq and 0.5 Mq stars using the PMS and 
main-sequence (MS) tracks as described in Section [2.21 In Fig. [10] the first few Myr are shown. Marked in the figure are also 
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logio(age/Myr) 

Figure 8. Logarithm of tiie number of stars versus logaritiim of tlic age for Model Fit 1. Solid histogram: The assigned input age. Dotted 
histogram: The recovered output age of all stars. 




-10 12 
logio(age/Myr) 

Figure 9. Like Fig. [S] but split for high- and low-mass stars. The dotted line shows the histogram of output ages for stars above 1 Mq 
and the long- dashed one for stars below that limit. 



the evolutionary tracks for four different single stars. A 1 Mq star, one of 1.25 Mq, sl 0.8 Mq and a 0.5 Mq star. Within its 
first 1.4 Myr the binary first looks like a 1.25 Mq star until it evolves parallel for some time to a 0.8 Mq star. The measured 
age of the binary is 0.11 Myr when it crosses the 1 Mq track, while the true age of the two companion stars at that moment 
is 0.05 Myr. Between 0.3 and 1.4 Myr the binary looks almost identical to a 0.8 Mq star but this star evolves through this 
part of its track at an age between 0.16 and 0.66 Myr. Therefore, the binary of 1.5 Mq, interpreted as a single star, looks for 
about 1 Myr like a 0.8 Mq star of only half its age. 

In Fig. [TT] the final contraction of the pre-main-sequence stars to the zero-age-main-sequence (ZAMS) and the first Myr 
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B-V 

Figure 10. The evolution of some young stars tlirougli tlie colour-magnitude diagram along the Hayashi contraction tracks for the first 
3 Myr. The dotted and short-dashed lines indicate the evolution of single stars, while the thick solid line shows a binary made of a 1 Mq 
and a 0.5 Mq star. The thin long-dashed lines represent isochrones for single stars with ages as indicated in the figure. For more details 
see the text in § 13.2.11 




B-V 

Figure 11. Like Fig. [TO]but following the PMS evolution until an age of 100 Myr. For details see the text in 13.2.11 



of the later MS evolution up to an age of 100 Myr is showijj. Here the binary evolves for some time parallel to the track of 
a 1.05 AIq star. But while the 1.05 Mq object moves through this phase when 9.5 to 22 Myr old the binary is about 5 Myr 
older, with a true age between 14 and 25 Myr. 



^ The 1.05 Mq star needs about 40 Myr to reach the ZAMS while the process takes roughly 80 Myr for the 1 Mq star. 
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3.3 Recovered total stellar mass 

As not all stars are recovered when observing unresolved multiple stars because they are either absorbed in pseudo stars or 
appear with a pseudo age too different to the bulk age, not only the ages are effected but also the deduced amount of mass 
in stars in a cluster. The two last columns in Tab. |3] show the fraction of recovered mass in the different models in proportion 
to the initial mass of the 10® original stars. First, all the "pseudo" masses for all recovered stars are added up and than 
divided by the sum of the real masses of all stars in the calculation. And second, only the "pseudo" masses of the stars with 
pseudo ages recovered within 25% of the real age are added up and compared with the real total mass of all stars. This second 
comparison is made because in a real observed sample stars with "wrong" ages might not be counted as cluster members but 
as field stars appearing superposed on the cluster. As can be seen in the second last column between 15 and nearly 60% of 
the true mass might be missed due to unresolved multiples. If only the stars are considered which are well within the true 
age of the systems (last column of Tab.|3J between 41 and 95% of the actual mass of the star cluster can be missed. It should 
be noted here that, as pointed out in § [1] the dynamics of massive stars in star clusters might change the mass- ratios of 
massive multiple systems towards more massive companions. This would result in an increasing amount of "hidden" mass in 
the cluster even if all stars were randomly paired in the beginning. This can be seen in Tab. [3] by the fact that both methods, 
"SpP" and "SyP", have lower recovered total masses than the "RP" method. Additionally, t he unresolved multiple stars ca n 
influence dynamical mass estimates of star clusters by biasing radial velocity measurements ijKouwenhoven fc de GriislboOSl l. 



4 DISCUSSION AND CONCLUSIONS 

A large numerical study has been carried out to determine whether unresolved multiple stars might hide a different underlying 
massive-star IMF in the observational data. In order to do so three different pairing mechanisms have been considered on 
order to produce the unresolved systems. 

The main difference between these three pairing mechanisms can be summarised as follows: Random Pairing (RP) keeps 
the MF of all stars and the MF of the primaries constant. As the MF of secondaries is constrained by the input IMF RP 
changes the system IMF. The "Special Pairing" (SpP) method with its mass-dependent lower mass limit for companions of 
massive stars keeps all MFs (above IMq) parallel but changes the system MF in a way similar to RP. As "System Pairing" 
(SyP) keeps the system MF as a constant input and puts constrains on the mass-ratio for massive systems, the primary and 
companion MFs are changed to fulfil these constrains. 

A first set of models (§ 13. 1|) has been studied under the assumption that all the masses of the stars and systems can be 
recovered exactly from the observations. While in the second set (§ I3.2|l the stars are evolved according to stellar evolution 
models, then merged following the different pairing algorithms and retrieved as if they were single stars. The results of this 
study of the effect of large numbers of unresolved multiple stars can be summarised as follows: 

• For random pairing and special pairing the system IMF is always slightly steeper above 1 Mq than the IMF for all stars 
(Tab. [21). 

• This, however, does not automatically lead to a steeper observed mass function. 

• In the case of random sampling there is n o effect on the observe d IMF above 1 Mq whatsoever (Tab. [3]). This seems 
to be in contrast with som e of the results by Mai'z Apellaniz ( 20081 ) . For certain combinations of lower and upper mass 



limits Maiz Apellaniz (|2008 ^) find differences in the slopes of up to Aa = ± 0.2. But when using a similar upper mass limit 
(120 Mq) as in this study (150 Mq) he finds only very little changes in the slope. He concludes that "in most cases the 
existence of unresolved binaries has only a small effect on the massive-star IMF slope". His conclusions are different for 
chance superpositions of stars which are not considered in this contribution. 

• For special pairing and system pairing the observed IMF is slightly steeper (Aa « 0.1 dex) than the IMF for all stars 
(see Fig. [12}. Therefore, an in reality observed IMF slope of 2.35 therefore might have a true underlying IMF for all stars of 
2.25 above a few Mq. 

• The effect is in general smaller than the error bars on observational slopes of IMFs. 

• Nonetheless, unresolved multiple systems can hide between 15% and up to 95% of the mass of a star cluster, depending 
on the pairing mechanism and the amount of higher order unresolved systems. 

• Furthermore, age dating in star clusters can by affected severely: 

- Even a single-age population will appear to have an age spread if the fraction of unresolved multiples is large. 

- Lower mass stars can appear 0.6 Myr older than massive stars in this case. 

Based on the broad range of models presented here we conclude that an observed Salpeter/Massey IMF slope (aobs = 
2.35) can only mask a slightly flatter IMF for all stars with atmc ~ 2.25. 
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Figure 12. Change in slope (Aa = Oobscrvcd ~ Cfallstars ) between tlie observed slope and a slope for all stars of the models with stellar 
evolution (see Tab. [3]l. Filled black dots indicate models with 100% binaries, open circles are the models with 100% triples and the 
"x'-symbols show models with 100% quadruples. Note that for Model Fit 4 (q = 1) only the 100% binary case is shown and that for 
Model Fit 2 (SpP) the values for the binaries and the quadruples are too close together to be distinguishable. 



It should be noted here that Mai'z ApeUaniz ( 20081 ) too finds little evidence for effects of unresolved, random paired, 
binaries on the slope of the IMF of massive stars despite using a different colour {U — V instead of B — V), different 
stellar evolution models, different routines to calculate the colours and another method to measure the slopes of the mass 
functions. But he further concludes that chance superpositions of stars in star clusters also introduce a bias on the slope 
determinati on, and th at the slope-fitting method via histogram construction additionally biases slope determinations (see also 
[Mai 'z Ape Uaniz fc tJbed a 2005). Due to our novel approach by using a bias- free maximum likelihood method (see § [2]3} to 
m easure the slopes of the mass fun ction s we do not suffer from the disadvantages of the histogram fitting method as described 
in lMai'z ApeUaniz fc Ubedal (|2005l ') and lMai'z ApeUaniz! \200i ). 
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APPENDIX A: THE STELLAR INITIAL MASS FUNCTION 



The following multi-component power-law IMF is used throughout the paper: 



, mil ^ ffi < tuq, 
, mo i^m < mi , 
,mi <m < rriniax, 



(Al) 



with exponents 

ao = +0.30 
ai = +1.30 
«2 = +2.35 
03 = +2.35 



0.01 ^ m/M© < 0.08, 
0.08 ^ m/Mg < 0.50, 
0.50 ^ m/MQ < 1.00, 
1.00 < m/MQ. < mma 



(A2) 



where dN = ^ (m) dm is the number o f stars in the mass interval m to m + dm. The exponents Qi represent the standard 
or canonical IMF ( Kroupa 2001 . 20021 ') . The advantage of such a multi-part power- law description is the easy integrability 
and, more importantly, that different parts of the IMF can be changed readily without affecting other parts. Note that this 
form is a two-part power-law in the stellar regime, and that brown dwarfs contribute about 4 per cent by mass only and 
need to be treated as a separate population su ch that the IMF has a discontinuity near mn ~ 0.08 AIq with k 



( Kroupa et al.ll2003 



Thies & 



Kroupal 2007 . 2008h . A log-normal form below 1 Mq with a power-law extension to h igh masses 



was suggested by Chabrton 200^ rbut is indistinguishable from the canonical form ( Dabringhausen et al. 20081 ) and does 



not cater for the discontinui ty. The observed IMF is today understood to be an invariant Salpeter/Massey power-law slope 
I Salpete il ll955l : lMassevll2003l ) above 0.5 Mq, being independent of the cluster density and metallicity for metallicities Z ^ 0.002 
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Figure Bl. The minimum possible secondary mass (^boundary) in dependence of rriprimaryi both in logarithmic units (cg. lBll l. Indicated 
by a dotted line is the lower limit for the secondary mass for a primary of 150 Mq. 



jMassev fc Hunterll 19981: 



Sirianni et alj 



200C, 



2OO3I : IPiskunov et al.ll2004l : IPflamm-Altenburg fc Kroupalbood ) 



2002!: barker et al.ll200ll : lMassevl[l99i . |2002| . booj IWvse et al.ll2002l : iBell et al 



The basic assumption underlying our approach is the notion tha t all stars in every cluster are drawn from this same 
universal parent IMF, which is consistent with observational evidence ( Elmegreen 19991 : Kroupa 2001). 



It should be noted he re that, while not indicated in eq. lA2l there is evidence of a maximal mass fo r stars (m 



IT 



150 M( 



Weidner fc Kroupall2004l ). a result confirmed by several independent studies IjOev fc Cl arke 2005; Figor 2005, 
and that a max;imal-star-mass-star-cluster-mass relation, minax(A^cci) ^■max*, exists jWeidner fc Kroupa, 2006 ) 



Koen 2OO1 



APPENDIX B: SPECIAL TREATMENTS TO PRODUCE SECONDARY STARS 
Bl Special Pairing (SpP) 

The following description has been developed in order to produce multiple stars for which the mass-ratio increases with 
increasing p rimary mass. Brown dwarfs (m < 0.08 M©) are not included as they are most likely to be treated as a separate 
population jThies fc KroupalboO?! . boosi ') . 



In the regime above 2 Mq the following quadratic relation was chosen to define the range of mass ratios as a fraction of 
the primary star mass, mprimary: the least massive allowed companion, mboundary, 

niboundary = 3.1 • 10"^ X ml^^^^^y + 0.034 X mprimary (Bl) 

Companions are paired randomly within the so defined mass range. This eq. assures that primaries of up to 2 Mq have a 
least-massive companion of 0.08 Mq, while the least-massive companion is 75 Mq for the most massive stars considered here 
(150 Mq). Fig. [BTI shows the minimum possible secondary mass (mboundary) in dependence of mprimary resulting from eq. IBll 
Stellar masses are created by randomly drawing all stars from the IMF. Binaries are then made from this set of stars as 
described in § 12.11 

The combined range of possible g- values for eq. IBll and random pairing for stars below 2 AIq is shown in Fig. IB2I The 
solid line marks the minimum allowed mass ratio. The companion masses lie within the shaded region. 

Fig. IB3I shows the resulting g-distribution for 4 different mass bins. While the lowest mass bin (0.08 - 2 Mq) shows a 
flat distribution resulting from random pairing the other three bins are peaked towards the lower edge of allowed mass-ratios. 
This is due to the steep power-law slope of the canonical IMF, which makes it highly likely that the chosen companion for a 
massive star is close to the lower limit mboundary. 
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Figure B2. The possible g- values {shaded region) in dependence of the primary mass (in logarithmic units) produced by the SpP 
procedure (Appendix IBH . Inside the shaded region the companions are randomly drawn from the IMF. The dashed line shows the 
minimum possible g- value for random sampling over all masses with a companion minimum mass of 0.08 Mq. 
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Figure B3. The mass-ratio (q) distributions for primaries in several mass intervals for the SpP model. Solid line: rviprimary from 0.08 
to 2 Mq, dotted line: 2-10 Mq, short-dashed line: 10 - 50 Mq and long-dashed line: 50 - 100 Mq. 



B2 System Pairing (SyP) 

In order to obtain systems with components which are strongly biased towards equal mass systems, the SyP scheme is moti- 
vated by the possible merger origin of massive stars and may also by applicable for the competitive accretion scenario. System 
masses are drawn randomly from the adopted fixed system IMF and then split into companions as follows: 
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Figure B4. The possible 5- values {shaded region) in dependence of the system mass (in logarithmic units) produced by the SyP procedure 
(Appendix IB2I I. But within the shaded region the q- values are only randomly distributed for msystcm < IOA/q. Above this mass the 
q- values are from distributions as shown in Fig. IB6I 

For msyatom > WMq: 

^^systcm X R 

' = 10 ' 

if q > 0.9 -.q = (0.1 x R) + 0.9; 

^^systcm 
JTlprimary = — ; ; ; 

1 + g 

'^companion — ^72systcm ^^primaryi 

where 7? is a random number between and 1. This formulation biases the mass ratio distribution for massive stars strongly 
towards more massive companions while still allowing for occasional low-mass secondaries. Systems less-massive than 10 Mq 
are split randomly in two or more halves with the limitation that the lower limit for the stellar mass of the companions is 
0.08 Mq. 

The IMF for the systems is the same as in §|A]but with different mass intervals. The lowest system mass allowed is 0.16 
Mq and the highest 300 Mq and the exponent changes at 1 Mq instead of 0.5 Mq. But for the stars produced from these 
systems the lower and upper mass limits are still 0.08 and 150 Mq. 

The combined range of possible g- values is shown in Fig. IB4I while the probabilities for the companion masses for a few 
examples are shown in Fig. IB5I These probabilities are chosen randomly from a Gaussian distribution between mboundary 
and the maximum possible mass for the secondary, msec, max, which is HisEiaa, The mean of the Gaussian is msec, max and its 
dispersion is given by a /(0.064msystcm +3.29). 

In Fig. IB6I the g-distribution is shown for 4 different mass bins. The three bins with stars above 10 Mq peak strongly 
above q — 0.8, indicating the strong tendency to produce nearly equal mass binaries with this algorithm. 




(B2) 
(B3) 
(B4) 

(B5) 
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Figure B6. The mass-ratio (g) distributions for primaries in several mass intervals (Appendix IB2I I. Solid line: mprimary from 0.08 to 2 
Mq, dotted line: 2-10 Mq, short-dashed line: 10 - 50 Mq and long-dashed line: 50 - 100 Mq. The mass-ratio distribution for low-mass 
stars {solid line) is not completely flat due to the lower mass limit for stars of 0.08 Mq. 
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Figure CI. IMFs of Models 1 to 4. The thick solid lines arc the IMFs of all stars, the thin solid lines arc the input IMFs (identical to 
the IMF of all stars for models 1 to 6), the thick dotted lines the system IMFs and the thin lines represent the mass functions for the 
individual components. The thin long-dashed lines are the primary star IMFs, the thin dash-dotted lines are the secondary star IMFs, the 
thin short-dashed lines are the tertiary star IMFs (only panels 2, 3, 5, 6, 8 and 9) and the thin long- dash- dotted lines arc the quartiary 
star IMFs (only panels 3, 6 and 9). The "spikes" in the system IMFs (thick dotted lines) in panels 4 to 10 are a result of changing the 
q-distribution at 2 Mq . 



APPENDIX C: IMFS WITHOUT STELLAR EVOLUTION 
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Figure Dl. IMFs of Models Fit 1 to Fit 3 from Table |3] Thick solid lines are tlie input IMFs, the thick dotted lines are the system 
IMFs without stellar evolution, the thick short-dashed lines are the "observed" IMFs, the long-dashed lines are the primary star IMFs, 
the dash-dotted lines are the secondary star IMFs, the thin short-dashed lines are the tertiary star IMFs (only Models Fit 2, 3, 5, 6, 8 
and 9) and the thin long- dash- dotted lines are the quartiary star IMFs (only Models Fit 3, 6 and 9). For the first panel ("Fl RP") the 
"observed" IMF is created only from recovered stars with pseudo ages less than 5 Myr. For each of the other three panels all recovered 
stars regardless of their ages are used. These panels are marked by "all ages". Note the dip in all models around logj^Q(m) -0.8 (about 
0.16 Mq). It is the result of differences in the PMS tracks from different authors used here. These differences shift some of the stars 
between 0.15 and 0.3 Mq to lower masses but have no influence on the results of this study as we restrain ourselves only to the high-mass 
end of the IMF. 



APPENDIX D: IMFS WITH STELLAR EVOLUTION 



This paper has been typeset from a T[h;X/ M?[5X file prepared by the author. 
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Figure D3. Like Fig.[DT]but for the Models Fit 7 to Fit 9 from Table|3] 
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